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Receptors of bacterial chemotaxis form clusters at the cell poles, where clusters act as "antennas" 
to amplify small changes in ligand concentration. Interestingly, chemoreceptors cluster at multiple 
length scales. At the smallest scale, receptors form dimers, which assemble into stable timers of 
dimers. At a large scale, trimers form large polar clusters composed of thousands of receptors. 
Although much is known about the signaling properties emerging from receptor clusters, it is 
unknown how receptors localize at the cell poles and what the cluster-size determining factors 
are. Here, we present a model of polar receptor clustering based on coupled trimers of dimers, 
where cluster size is determined as a minimum of the cluster-membrane free energy. This energy 
has contributions from the cluster-membrane elastic energy, penalizing large clusters due to their 
high intrinsic curvature, and receptor-receptor coupling favoring large clusters. We find that the 
reduced cluster-membrane curvature mismatch at the curved cell poles leads to large and robust 
polar clusters in line with experimental observation, while lateral clusters are efficiently suppressed. 
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I. INTRODUCTION 

Chemoreceptor clustering is widely conserved among 
bacteria and archaea allowing cells to detect chemi- 
cals in the environment with high sensitivity over a wide 
range of background concentrations. In the bacteria Es- 
cherichia coli, Salmonella enterica, and Caulobacter cres- 
centus, receptor clustering is well documented and occurs 
at multiple length scales. At a small scale, chemotaxis 
receptors form stable homodimcrs, which then assemble 
into larger complexes with receptors of different chemi- 
cal specificities intermixed Q. Three homodimers, con- 
nected at their signaling tip, form a trimcr of dimers 




believed to be the 
At a larger scale, 



(named trimer from here on) 
smallest stable signaling unit 
thousands of receptors [7[ form approximately 200-nm 
large polar clusters (cf. Fig. la) |[ QJ, El EI EE El • 
Despite the excellent characterization of much of the bac- 
terial chemotaxis network, it is unknown how receptors 
localize at the cell poles and how they assemble into large 
polar clusters. 

Polar localization appears to be an intrinsic property 
of chemoreceptors [UEII ■ It hardly depends on the pres- 
ence or absence of the receptor-bound kinase CheA and 
adapter protein CheW [l|| , and is unaffected by removal 
of the periplasmic ligand-binding domain of the receptors 
[l9| . It is also a passive process since newly synthesized 
receptors, initially inserted at random positions in the 
membrane, diffuse and ultimately become trapped at the 
cell poles [IE Hi) . Most importantly, polar localization 
appears to depend on membrane curvature. First, in- 
hibition of actin-homologuc MreB in growing cells leads 
to cell swelling and a diffuse receptor distribution, with 
remaining receptor localization in areas of increased cell 
curvature p2| . Second, receptor-membrane extracts self- 
assemble into round micelles after receptor overexpres- 



sion and cell lysis [23[ (Fig. lb). From electron micro- 
graphs, the intrinsic curvature of the trimer structure can 
be estimated 2J]. Third, other two-component receptor 
dimers, e.g. the receptor LuxQ of the quorum-sensing 
pathway, dimerize without forming trimers of dimers and 
are evenly distributed over the cell surface 25j . Taken to- 
gether, these observations suggest that the distinct trimcr 
structure with its increased intrinsic curvature is respon- 
sible for polar receptor localization. 

While trimers may have a tendency to localize at 
the cell poles and areas of high membrane curvature, 
tight clustering requires an attractive coupling among 
the trimers. The conventional view is that CheA and 
CheW mediate interactions among receptors. Alterna- 
tive models include swapping of the cytoplasmic recep- 
tor domains [2(| and membrane-mediated coupling (2?| 
(see Results and Discussion section). The high sensitiv- 
ity and cooperativity obtained from in vivo FRET (fluo- 
rescence resonance energy transfer) l29j and in vitro 
[30T ] data demonstrate that the functional units of recep- 
tor signaling are indeed larger than trimers. These ob- 
servations are suppor ted by recent quantitative models 

m, mL us, m, hi, in . 

Based on these observations, we propose a model for 
polar receptor localization and clustering due to the high 
intrinsic curvature of trimers and an attractive trimer- 
trimer coupling. Specifically, we consider a membrane- 
embedded cluster composed of trimers. For a sphero- 
cylindrical cell, we assume that the average membrane 
curvature at the poles is twice as large as average cur- 
vature at the lateral surface area, and that trimers have 
a high intrinsic curvature (Fig. 2a). The intrinsic cur- 
vature of a trimer tends to deform the membrane (Fig. 
2c), penalizing large clusters of trimers. However, at- 
tractive coupling between trimers favors cluster forma- 
tion, leading to a competition between these two op- 
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FIG. 1: Electron micrographs of chemoreceptor clusters, (a) 
Extended clusters in membrane preparations [1, 0, HcL [lj. 
Clusters of similar size are observed at the poles of living 
cells P2, 0- Image is courtesy of Michael Manson. (b ) 
Self-assembled round micelle at receptor-dimer resolution [23| . 
Image printed with written permission from publisher. All ex- 
periments are based on Tsr-receptor overexpression, as well 
as membrane extraction, negative staining, and freezing for 
imaging. 



posing energy contributions. Using continuum elastic 
theory, we derive an analytical expression for the to- 
tal cluster-membrane energy We find that, due to the 
reduced cluster-membrane curvature mismatch at the 
poles, trimers favorably cluster at the poles and not at 
the lateral cell area. Furthermore, the cluster-size distri- 
bution is determined by the cluster-membrane energetics, 
as well as the trimer density in the cell membrane. Our 
predicted average cluster size is in line with experimental 
observation. 



II. MODEL 



Receptor geometry 

In our model, receptor dimers are assumed to always 
be associated in trimers, the smallest stable signaling 
unit 0, @- In Fig- 2a, the receptor dimer length and 
width, as well as the distance between neighboring 
dimers within a trimer are taken fro m p artial crystal 
structures [1] and electron microscopy (23[. Very similar 
parameters were used to model the physical response 
of trimers to osmolytes measured by homo-FRET [36j . 
Importantly, the estimated value for the intrinsic curva- 
ture of a trimer corresponds closely to the inverse radius 
of self-assembled micelles (Fig. lb, HH). The value for 
the trimer cross section A (table I) is consistent with a 
three-dimensional model of the receptor cluster 37l 38j 



and estimates from cryo-electron microscopy [1, [R3 . 123 1 



Elastic cluster-membrane energy 

The elastic energy of a membrane-embedded cluster is 
determined by cluster and membrane bending energies 



a Trimer of dimers 
10 nm 



b Cell pole 
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inner 
membrane 



turgor pressure 
C Lateral cell area 



FIG. 2: Schematic of membrane-inserted receptors. (a) 
Trimer (orange shaded area) of dimers (red bars) with geomet- 
ric parameters. Ct is the intrinsic curvature of a trimer. (b) 
Cluster of three trimers at the cell pole. Trimer-trimer cou- 
pling strength Jo is indicated by green arrows. Also shown are 
the cell wall and the inner membrane (gray) with curvature 



The height profile h{r) describes the cluster-membrane 



deformation as measured relative to the preferred height due 
to cell wall and turgor pressure, (c) Same cluster at a lateral 
position with membrane curvature C m ,i- Red arrows indicate 
energetically unfavorable membrane deformations due to the 
cluster-membrane curvature mismatch. 



and a pinning potential 



E,. 



^(2C c (f) - C T f + ^(h c (r) - h (f)) 2 



d z r 



^(2C m (r) - C m f + ^(h m (f) - h a (f)) 2 



d 2 r. (1) 



The first term, proportional to the bending stiffness k c 
of the receptor cluster, penalizes deviations between the 
total cluster curvature 2C C and the preferred cluster cur- 
vature, which is equal to the intrinsic trimer curvature 
Ct- The total cluster curvature 2C C = C1+C2 is defined 
by the two principal curvatures C\ and Ci [3!| . The sec- 
ond term in Eq. [TJ proportional to the pinning modulus 
A [4(| El| , penalizes deviations of the cluster height h c 
from the preferred height ho determined by the shape of 
the curved cell wall. The third and forth terms in Eq. [1] 
mirror the first and second, respectively, and describe the 
cluster-surrounding membrane with total curvature 2C m , 
preferred curvature C m , and height h m . Bending stiff- 
nesses k c and from optimal packing of receptors 
and lipids, respectively, aiming to protect hydrophobic 
residues from polar water. The pinning modulus arises 
due to the turgor pressure, which pushes the membrane 
and cluster outward, while the rigid cell wall pushes them 
inward. The net effect is a penalty for deformations away 
from the preferred cell shape (see Fig. 2) [4fj|.[4l|. 
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Let us define the relative height perturbation h c r m ) = 
h c ( m ){r) — ho(r) as measured relative to the preferred 
height h[)(f) in the direction of the normal, pointing in 
the radial direction outward from surface ho(f) (Fig. 
2 b,c). This allows us to perform the following cal- 
culations in the so-called normal gauge, where the to- 
tal curvatures for the cluster (and membrane) 2C c ( m ) 
arc now given by C m + (C m - 2C G )h c(m) + V 2 /i c(m) to 
lowest order, following from the first-order variation of 
the geometry j^]. Here, Cq is the Gaussian curvature. 



We use two further approximations: (1) V 2 /i c ( m ) :» 
(C 2 ^ — 2CG)/i c ( m ), justified for small amplitude and large 
bending ripples. Hence, the contribution proportional 
to /i c ( m ) is neglected. (2) The curvilinear Laplacian is 
replaced by the flat Laplacian V 2 « d 2 /dx 2 + d 2 /dy 2 , 
valid for sufficiently small clusters. Introducing 2C c ( m ) ~ 
C m + V 2 /i c ( m ) and AC = CV — C m . Eq. [T]can thus be 
written as 



J 



E, 



^ { V 2 h c (r)-AC) 2 + ±h 2 c (r) 



d 2 r + 


I 




J m 



r 



^ 2 h m (r)) 2 + \h 2 m {r) 



d 2 r. 



(2) 



Parameter AC is the important curvature mismatch, i.e. 
the difference between the cluster and membrane curva- 
tures. The elastic energy model in Eq. [2] neglects sur- 
face tension, as well as an Gaussian curvature effect due 
to the different cluster and membrane elastic properties 
(see Appendix A for a justification). In this paper, the 
model is applied to both polar and lateral clusters using 
standard parameter values given in table U Very similar 
models were previously applied to describe mixtures of 
lipids with different curvatures [H, EH, EE 0, EH ■ For a 
general review and alternative elastic energy models, see 

El, Ell. 

Considering a circular cluster of radius R, the total 
height profile h(r) = h(x,y) is composed of h c (f) for the 
cluster (r < R) and h m (f) for the membrane (r > R), 
and is determined by minimizing the total elastic cluster- 
membrane energy with respect to variation of h(x,y). 
Following Ref. |48[, minimizing an clastic energy of the 
generic form 



b rd T f , dh dh d 2 h d 2 h\ 

* [x,y,h, — , — )dxdy (3) 



dx ' dy ' dx 2 ' dy 2 

leads to the following Euler-Lagrangc equation 

9* d ( 9* \ 5/9* 
~dh ~ dx \d(dh/dx) J ~ dy~ \d(dh/dy) 
d 2 ( 9* \ d 2 ( 9* 



dx 2 \d(d 2 h/dx 2 



dy 2 \d(d 2 h/dy 2 ) 



(4) 
0.(5) 



Replacing * by the expressions in Eq. [2] leads to the 
forth-order linear differential equation [i^, [13, HH for 
cluster (membrane) 



V 4 h, 



(m) 



A 



c(m) 



■h c ( m ){r) = 



(6) 



which is independent of AC in the small deformation 
approximation. Since we consider a circular cluster lo- 



cated at the origin of the coordinate system, we ap- 
ply cylindrical symmetry from now on. To solve Eq. 
[5] for cluster and membrane one needs four boundary 
conditions for each equation. We require dh c /dr = 
at the origin and that the membrane deformation van- 
ishes far away from the cluster, i.e. lmv^oo h m (r) = 
and linv^oo dh m (r)/dr = 0. We further impose that 
the solutions for cluster and membrane match at the 
cluster-membrane interface, i.e. h c (R) = h m (R) and 
dh c /dr\ R = dh m /dr\ R . 

Eq. [5] can be solved by applying the Kelvin differential 
equation 



S7 2 h(Pr) - i(3 2 h{(3r) = 



(7) 



leading to /3 c ( r 



for the cluster (membrane). 



In Eq. the second-order derivative is now calculated 
using V 2 = d 2 /dr 2 + 1/rd/dr. The solution to Eq. [7] is 
given by [52| 

ber ((3 c(m) r) + ibei (f3 c{m) r) = I {(3 c{rn) re- i37r/4 )(8) 
ker (/3 c{m) r) + ikei (/3 c{m) r) = K (f3 c{m) re l7T/4 ) (9) 

where berg, beio, fcero, and keio are the zeroth-order 
Kelvin function, and Jo and Kq are the zeroth-order mod- 
ified Besscl functions of the first and second kind, respec- 
tively. 

To construct the solution for the cluster (r < R), only 
the modified Besscl function of the first kind has zero 
slope at r = 0, whereas to construct the solution for 
the membrane (r > R) only the modified Bessel function 
of the second kind has a vanishing real part at infinity. 
To obtain real solutions, we need to add their complex 
conjugates 

h c (r) = {a + ib)I {f3re- l3 ^ 4 ) + (a - ib)I Q {f3re + ^/ 4 llO) 
h m ir) = (c + id)K Q ((3re+' l ^ 4 ) + (c - id) K ((3r e^^ll) 
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where a, 6, c, and d are real parameters to be determined 
by matching the boundary conditions. 

Solution coefficients 

After matching boundary conditions, the coefficients of 
the solution in Eqs. [TOl and [TT1 are given by 



Ho 
2 

23Jo(/3 c ,2-R) 
Hi + 2d%K (p mA R) 

S$tK (J3 m ,iR) + H 1 St\p m ,iK 1 (J3 m ,iR)] 
2Qf\p m ,iK (p mi i)K 1 (p mil R)] 



(12) 
(13) 
(14) 
(15) 



where S = 25ft[(a + ib)f} c ^h{(3 c ,2 R )\ is the slo P e 
at the cluster-membrane interface, Hq is the clus- 
ter height at r = 0, and Hi is the height at 
r = R, i.e. at the cluster-membrane interface. 
We further used dIo(/3r)/d(f3r) — f3I\((3r) and 
dKo(f3r)/d((3r) = — /3_Ki(/3r). Remaining unknown 
parameters, such as Hq and Hi, are determined by 
numerically minimizing the elastic energy. 



Analytic expression for elastic energy 

The integrals in Eq. [2] can be solved analytically using 
integration-by-parts twice 0, [53[ and exploiting Eq. [6] 
leading to 



E el = ttk c R{S ■ V 2 h c \ R - Hi ■ V 3 h c \ R - 2AC ■ S} + -n Kc AC 2 R 2 - nK m R{S ■ V 2 h m \ R - Hi ■ V 3 h m \ R }. (16) 



r 



The higher-order derivatives are calculated from the so- 
lution Eqs. ITD1 and [TT1 using [52j 



VJ (/?r) 
V 2 / (/3r) 

V 3 / (/?r) 



— + -—) I (/3r) 
dr 2 r dr J 



(17) 



(18) 



dr 



V 2 / (/3r) 



dr 3 r dr 2 r 2 dr 



and 



VK (l3r) = -pKiffir) 
V 2 K (Pr) - [3 2 K {l3r) 
V 3 K {f3r) = -p 3 Ki(pr). 



(19) 



(20) 
(21) 
(22) 



In the limit of very large clusters, the elastic energy 
Eq. [TrJl reduces to 



E 



(oo) 



^ttk c AC 2 R 2 = ^k c AC 2 AN, (23) 



since the other contributions to the elastic energy 
grow more slowly with size (c/. Fig. 5a). The number 
of trimers in a cluster of radius R is given by N w R 2 ir/A. 



energy of a cluster of radius R is described by 

E & = -J(N) ■ N. (24) 

For a triangular lattice used in Fig. 3, we consider the 
following two expressions for the average coupling energy 
per trimer. Assuming a cluster made of concentric rings 
formed around a central trimer, the coupling energy is 
given by [4l| 



(25) 



where Jo is the coupling energy between two neighboring 
trimers (table I). This expression is exact in the limit of 
large clusters. Alternatively, we fit 



J(N) = 3 Jo 



N 



N + N Q . 5 



(26) 



to exact interaction energies of small compact clusters, 
and obtain parameter iVo.5 = 4.72. According to Fig. 
3, Eq. overestimates the interaction energy for small 
clusters, whereas Eq. [26] overestimates the interaction 
energy for large clusters. In the limit N — > 00, both 



models produce E : 



(00) 



-3JoiV. In the following we 



use Eq. [26] since it helps stabilizing large clusters. 

The total energy is the sum of elastic energy and the 
attractive energy 



E(N) = E el {N) + E a (N), 



(27) 



Attractive trimer-trimer coupling 

In our model, trimers interact favorably when in close 
contact, driving cluster formation. The total coupling 



where we explicitly included the dependence on the 
cluster size N. 
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Distribution of cluster sizes 

For a sphero-cylindrical cell, the preferred membrane cur- 
vature at the poles C m — C m . p is twice as large as the 
preferred membrane curvature at the lateral area C m — 
C m ,i = Cm,p/2, however smaller than the trimer curva- 
ture Ct- This leads to a smaller cluster- membrane cur- 
vature mismatch at the poles AC P = Ct — C m ^ p than at 
the lateral area AC/ = Ct — C m> i = Cy— C m>p /2. Conse- 
quently, the total energy Eq. [57]of a cluster is also smaller 
at the poles than at the lateral area E P (N) < Ei(N), fa- 
voring polar clustering. 

Based on the total energies, statistical mechanics is 
used to calculate the cluster-size distribution at the 
poles/lateral area (p/l) [ill. 



P p/l (N) = Ne-^/iW-N^/k, 



(28) 



where the exponential Boltzmann factor describes the 
probability to observe a cluster of N trimcrs at the 
poles/lateral area for chemical potential /i. The chem- 
ical potential represents the energy required or released 
by inserting a trimer into the membrane, and is adjusted 
to fulfill an overall target trimer density on the cell sur- 
face (occupancy fraction) p via 



J2[P P (N) + Pi(N)] = p. 



(29) 



TABLE I: Summary of standard parameters. These 
parameters are used throughout calculations unless specified 
otherwise. 
Parameter Value Meaning 



A [tfiT/nm 4 
n c [k B T] 
n m [k B T] 
C~\ [nm 1 
C 



-l 

If 

Crp 

J [k B T 



nmj 
nml 



0.25 [40] Pinning modulus 
120 * Bending stiffness of cluster 
25 [40] Bending stiffness of membrane 
400 [39] Inverse of polar membrane curvature 
800 ' Inverse of lateral membrane curvature 
37 [231 Inverse of trimer-of-dimer curvature 
3 * Trimer-trimer coupling strength 
A [nm 2 ] 200 [37] Trimer cross section 
* these parameters are varied in Fig. 5C to check for 
robustness. ' based on a sphero-cylindrical cell. 

poles and lateral area is provided by 



Ae = 



(oo) (oo) 

—^-(Ct — 3/4 C m ^ p )C m , p , 



(32) 



N 



where we used C„ i; / = C„ hp /2. Increasing Ct beyond 
3/4C mjP favors polar over lateral clusters. Specifically, 
for iVAe > lk B T, a cluster of N trimers is significantly 
more favorable at the poles than at the lateral area at 
temperature T. 



Using the cluster-size distributions, the average cluster 
sizes at the poles and lateral area are given by 



(N) p/1 =J2 NP p/i( N )- 



(30) 



N 



Conditions for large polar clusters 

To find the conditions which favor polar clustering, we 
consider the total energy density e = E p /i/N, i.e. the to- 
tal energy of the membrane-embedded cluster per trimer. 
Generally, minimization with respect to N determines 
the energetically preferred cluster size. We note the fol- 
lowing. First, the energy density is generally a mono- 
tonically decreasing function of N, which eventually sat- 
urates for large N (cf. Fig. 5a). This indicates that 
maximal cluster sizes are energetically favorable. Sec- 
ond, the elastic energy density is always smaller at the 
poles than at the lateral area, demonstrating that polar 
clustering is energetically predominant. 

Let us consider the total energy densities in the limit 
N — > oo. In this limit, the total cluster- membrane energy 
density at the poles/lateral area is given by 



- ( 2 KcAC p/ iA ~ 3J ° 



(31) 



Consequently, the energy-density difference between the 



III. RESULTS AND DISCUSSION 

Based on experimental observations outlined in the In- 
troduction section, we propose a model for polar recep- 
tor localization and clustering (Fig. 2). The ingredients 
and model assumptions are as follows: (1) An individ- 
ual trimer of dimers (trimer), believed to be the smallest 
stable signaling unit 0, @, has a high intrinsic curva- 
ture Ct (Fig. 2a). The cell membrane has a higher 
curvature at the cell poles C m . p than at the lateral area 
C m> ;. For a sphero-cylindrical cell, we have specifically 
C m ,i = Cm,p/2. Since C T > C m , P > C m ,j, individual 
trimers favor the cell poles energetically, although this 
effect is very small by itself (fraction of thermal energy 
k B T). (2) Trimers are coupled with strength Jo when in 
close proximity (Fig. 2bc), driving cluster formation at 
the poles and lateral area (Fig. 3). (3) Due to the cluster- 
membrane curvature mismatch, growing clusters deform 
the membrane and are energetically penalized. However, 
since the cluster-membrane curvature mismatch at the 
poles AC P = C't — C m ,p is smaller than the correspond- 
ing mismatch at the lateral area AC/ = Ct — C m>p /2, 
this energy penalty is reduced at the poles (Fig. 2b). As 
outlined in the Model section, the model is implemented 
by considering a membrane-embedded cluster of radius 
R. The height profile of the cluster and the membrane 
minimizes the elastic energy, which is determined by the 
cluster and membrane preferred curvatures (respective 
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Cluster size N 




10 15 20 25 
Distance r(nm) 
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FIG. 3: Average interaction energy J(N) per trimer (in units 
of coupling strength Jo) as a function of cluster size (number 
of trimers N). Trimers are shown by orange disks, trimer- 
trimer coupling of strength Jo is shown by green bars. Clus- 
ters are assumed to have a triangular lattice structure. Sym- 
bols: small compact clusters of minimal circumference. Solid 
line and dashed lines: "Ring model" and "Fit-to-symbols 
model", respectively (see Model section). Horizontal dotted 
line: average interaction energy for infinitely large cluster. 
Inset: Possible mechanisms of trimer-trimer coupling includ- 
ing coupling mediated by receptor-bound CheA-CheW and 
membrane deformations based on large hydrophobic trans- 
membrane domains (TMs). 



Ct and C m p /;) and their bending stiffnesses (respective 
k c and K m ). Furthermore, a pinning modulus A |40l. l4lj 
pushes the membrane and the cluster against the rigid 
cell wall (Fig. 2 be) . This penalizes large deformations of 
the cluster and the membrane. The main findings are as 
follows: (1) Considered separately, poles and lateral area 
favor maximal cluster sizes energetically. (2) Actual clus- 
ter size is determined by timer density (entropy), where 
increasing trimer density pushes distribution of cluster 
sizes to larger values. (3) Polar-only clustering is a result 
of the reduced curvature mismatch at the poles, energet- 
ically stabilizing polar clusters and suppressing lateral 
ones. Our results are in line with the experimental ob- 
servation of large polar clusters, which were found to be 
robust [IlJMr 



J19( and only slightly affected by attractant 
binding 55l" 56] . expression level variation [H, H3]i 
and receptor methylation [ll], [H, [H, [6(| . 

Fig. 4 shows typical cluster- membrane height profiles 
for two different cluster radii, i?=10 and 100 nms. The 
height profile minimizes the cluster-membrane elastic en- 
ergy given in Eq. O The profile of the small cluster 
bulges out into the periplasmic space in a convex man- 
ner, while the large cluster of physiological size is flat- 
tened. The latter effect appears consistent with images 
from cryo-electron microscopy [IH, [l3|. In our model, 
large deformations are suppressed by the pinning poten- 
tial (left inset in Fig. 4b). Significant reduction of the 
pinning modulus A leads to strongly curved clusters (right 




Distance x(nm) 



Distance y(nm) 



FIG. 4: Cluster-membrane height profiles, (a) Profile h(r) 
as a function of distance r from center for a small cluster of 
radius R = 10 nm. Inset: Profile for large cluster of radius 
R = 100 nm. (b) Three-dimensional profile for small cluster. 
Left inset: Same for large cluster. For parameters k c , K m , and 
A see table I. Right inset: large cluster for very small pinning 
modulus (A = 10" 5 k B T/nm 4 ). 



inset in Fig. 4b). Note that the maximal deformation of 
the large curved cluster can easily exceed the width of 
the periplasmic space (20 nm [6ll]). emphasizing the im- 
portance of the pinning modulus. 

How large are clusters and how do their sizes differ 
between cell poles and lateral positions? Consider a 
membrane-embedded cluster of radius R or number of 
trimers A" ps R 2 it/A, where A is the trimer cross section. 
The total trimer energy E is equal to the sum of the un- 
favorable elastic energy E e i and the favorable attractive 
energy E a . Dividing by the number of trimers A" results 
in the corresponding energy densities e, e e i, and e a . Gen- 
erally, the minimum of the total energy density e as a 
function of cluster size N provides the preferred cluster 
size. Fig. 5a shows that key requirements for stable po- 
lar clusters are fulfilled: (1) The energy density reaches 
its lowest value in the limit A" — > oo, energetically favor- 
ing maximal clusters. (2) Although the energy-density 
difference Ae between the poles and lateral cell area can 



7 



be smaller than the thermal energy fceT, a cluster of N 
trimers is stabilized at the poles and suppressed at the 
lateral area when NAe is larger than UbT (see Model 
section). 

Due to the finite trimer density, cluster sizes are al- 
ways finite. To obtain the predicted distribution of clus- 
ter sizes, we consider the combined system of cell poles 
and lateral area. The distribution of cluster sizes can 
be calculated using Boltzmann statistics of equilibrium 
statistical mechanics (see Model section). A chemical 
potential is further adjusted to obtain a certain target 
trimer density. Fig. 5b shows the size distributions of 
polar and lateral clusters for three different trimer den- 
sities. Even at low trimer densities, very few residual 
trimers, stabilized by entropy, remain unclustcrcd at the 
poles and lateral area. The average radii of polar clusters 
shown in Fig. 5c correspond well with the observed clus- 
ter diameters of about 200 nm, whereas lateral clusters 
are significantly suppressed. 

As shown in Fig. 5c, our model predicts the aver- 
age cluster size as function of trimer density. This pre- 
diction can experimentally be tested through imaging. 
The fluorescence intensity of a cluster, e.g. measured 
using receptor-GFP fusion proteins [2(|, is proportional 
to the number of receptors in the cluster. Alternatively, 
imaging by cryo-electron microscopy can provide spatial 
cluster dimensions fl2l . [l3j . The dependence on trimer 
density can be studied by expressing receptors from an 
inducible plasmid. An recent experiment using the fusion 
protein CheY-YFP as a fluorescent marker indeed indi- 
cated a strong correlation between r ecep tor expression 
level and polar fluorescence intensity 62j. Furthermore, 



the predicted link between membrane curvature and clus- 
tering can be tested by quantifying receptor-fluorescence 
intensities for cell-shape phenotypes, e.g. when inhibit- 
ing actin-homologue MreB [2(| , responsible for rod shape 
in bacteria. Alternatively, cocci cells or round membrane 
vesicles can be used, allowing the study of receptor clus- 
tering in presence of only a single membrane curvature. 
In Fig. 6, we show the predicted average cluster size as 
a function of coccus radius and receptor density (expres- 
sion level). Increasing the coccus radius decreases the 
coccus curvature and leads to a larger cluster-membrane 
curvature mismatch, which reduces cluster size. In con- 
trast, increasing the receptor density shifts cluster size 
distribution toward larger clusters. 

In recent experiments the physical response of dimcrs 
were measured by homo-FRET using receptor- YFP fu- 
sions [H, [63[. These data indicate that the dimer- 
dimer distance in a trimer (distance between receptor 
C-termini) shrinks by 10% upon osmolyte stimulation. 
Osmolytes act as repellents and are presumably sensed 
through receptor-membrane coupling. To see if polar 
clustering is robust against such perturbations, Fig. 5c 
shows that a 10% increase of the dimer-dimcr distance 



(C 7 



36.2 nm, lower green error bars) destabilizes 



clusters very little, wheres a 10% decrease of the dimer- 
dimer distance (C^ 1 = 38.8 nm, upper green error bars) 
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FIG. 5: Quantifying polar receptor clustering, (a) Cluster- 
membrane energy density, i.e. energy per trimer, as function 
of cluster size N. Total energy density e (black) is sum of re- 
pulsive elastic energy density e e i (blue) and attractive cluster 
energy density e a (red). Solid lines correspond to cell poles, 
dashed lines to lateral positions. Thin solid and dashed blue 
lines indicate asymptotic limit of the elastic energy density 
for N — > oo at the poles and lateral positions, respectively. 
Polar clusters are stabilized by energy density Ae. (b) Distri- 
bution of cluster sizes for three values of the trimer density 
(occupancy fraction) p = 0.4, 0.1, and 0.025 at the poles (solid 
lines) and lateral positions (dashed lines), (c) Average clus- 
ter radius {R) at the poles (black solid line) and lateral area 
(black dashed line) as a function of trimer density p, based 
on standard parameters, including trimer-trimer coupling Jo, 
trimer curvature CV, and cluster bending stiffness k c from 
table I. Error bars indicate robustness to parameter changes. 
Upper error bars for poles and lateral area: 1.1 Jo (red), 0.98 
Ct (dimer-dimer distance reduced by 1 nm, green), 0.92 k c 
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FIG. 6: Predicted receptor clustering in cocci cells. Average 
cluster radius (R) as a function of coccus radius. Solid lines of 
varying thickness correspond to the three densities from Fig. 
5b. 



stabilizes clusters even further. To illustrate the robust- 
ness of polar clustering with respect to other model pa- 
rameters, we also varied cluster bending stiffness k c and 
trimer-trimer coupling strength Jq. Fig. 5c shows that 
these parameter variations only lead to modest changes 
in cluster stability (blue and red error bars, respectively). 

While our model provides a robust mechanism for the 
formation of large polar clusters, there are limitations of 
the model. (1) A recent study showed that, in addition to 
large polar clusters, there are also small lateral clusters 
at future division sites, such as 1/2, as well as 1/4 and 
3/4 cell length Q. However, these lateral clusters ap- 
pear immobile, presumably due to anchoring, and hence 
may form through a different mechanism. In Rhodobac- 
ter sphaeroides immobile lateral clusters of chemotaxis 
receptor homologues were even found in the cytoplasm 
[64] ]. Our model does suggest cluster formation at the 
new poles once cell division occurs and the membrane 
pinches off. Newly synthesized receptors, inserte d by 
the Sec-machinery throughout the cell surface [13, 165} . 
would begin to cluster at the new cell poles. However, 
if equilibration is too slow to grow new polar clusters 
from scratch after cell division, lateral clusters may be 
useful by serving as nucleation sites. (2) Structural work 
on receptors and receptor-bound proteins in Thermotoga 
maritima suggests that, at least for this bacterium, re- 
ceptor dimers may assemble into linear oligomers and 
not trimer-based clusters [H^]. A bioinformatics study 
on chemotaxis receptors across many species also ad- 
dresses this issue, but does not favor one model over the 
other [67j . On the other hand, cryo-electron microscopy 
of Caulobacter crescentus strongly supports trimer-based 
clusters [l3}. (3) Our model does not include interactions 
between multiple clusters. Such interactions may largely 
be unimportant, since fluorescence images generally in- 
dicate only one, rarely more clusters per cell pole [li} . 
(4) It can further not be ruled out that receptors do not 



localize to the poles themselves, but that certain polar 
lipids, e.g. cardiolipin [68} , provide favorable sites for 
receptor localization and clustering. (5) Although some 
clastic properties are included in our model, others are 
not, e.g. membrane thickness deformation due to large 
hydrophobic transmembrane domains [49] . Thickness de- 
formation leads to a line tension, i.e. an elastic energy 
proportional to the cluster circumference 2?tR, affecting 
both polar, as well as lateral clusters equally. While this 
effect does not change the stability of polar clusters, the 
line tension may provide a mechanism for trimer-trimer 
coupling [r3^ | (see next paragraph) and have influence on 
the cluster size distribution. 

What are the possible mechanisms responsible for 
trimer-trimer coupling? (1) Coupling mediated by CheA 
and CheW (Fig. 3, inset). Presence of CheA and CheW 
increases polar clustering only modestly [HI, Qi| [59} . 
However, overcxprcssion of CheA decreases, whereas 
overexpression of CheW increases the receptor cooper- 
ativity measured by FRET [29}. (2) Coupling medi- 
ated by elastic membrane deformations. Receptor ac- 
tivity has been shown to depend on receptor-membrane 
interactions [H, [63|, [7(j , which, we speculate, may pro- 
vide a mechanism for trimer-trimer coupling [69i . 171} . In 
support of this mechanism, receptor transmembrane re- 
gions are unusually large (24 to 30 residues) compared 
to the membrane thickness (30A or about 20 residues) 
[72} . Such large regions possibly lead to significant mem- 
brane deformations to protect the hydrophobic receptor 
residues from water (Fig. 3, inset). If transmembrane re- 
gions changewith activity, e.g. within the receptor piston 
model [7fj| l73l. 1741 [75} , trimer-trimer coupling may even 
depend on the receptor activity state, as proposed for 
the approximate two-state osmolarity-sensing MscL pore 
[27} . (3) Coupling mediated by swapping of cytoplasmic 
domains of neighboring dimers (26} . 

Many other sen sory receptors cluster as well, including 
B-cell W, T-cell |z3, FC7 [z3, synaptic [zl, and ryan- 
odine [80} receptors. This indicates that receptor cluster- 
ing is an important regulatory mechanism of the cell, e.g. 
to adjust signaling properties, recruit auxiliary proteins, 
or kinetically proof-read unexpected stimuli. Unlike bac- 
teria, cukaryotic receptor clustering appears much more 
dynamic, including receptor diffusion and internalization, 
and the underlying physical mechanism remain little un- 
derstood. Our work adds additional support to the idea 
that the elastic properties of receptors and membrane 
may be a general design principle to regulate receptor 
localization and clustering [47], [Slj] . 



APPENDIX A: SURFACE TENSION AND 
GAUSSIAN CURVATURE 

The elastic energy in Eq. [3] neglects surface tension 
and Gaussian curvature terms. Here we show that these 
two elastic energy contributions are indeed very small. 
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Surface tension 

In the small deformation approximation (Monge repre- 
sentation) for cluster and membrane, the contribution 
from the surface tension to the elastic energy Eq. [5] is 
given by 



E. 



el = £ / (vh c yd<r+ ^ / (Ai) 



where <J c (m) 1S the surface tension of the cluster (mem- 
brane) and V = (d/dx,d/dy). Surface tension arises 
in part from the attractive receptor-receptor and lipid- 
lipid interactions. Minimization of the total elastic en- 
ergy (Eq. [2] plus Eq. IA1[) leads to the Euler-Lagrange 
equation 

V% (m) (r) - ^±V 2 h c[m) (r) + _*_fe c(m) (f) = 

^c(m) ^c(m) 

(A2) 

for the cluster (membrane). Using integration- by-parts 
twice [49], [53[ and Eq. IA2i the elastic energy contribution 
due to the surface tensions of the cluster and membrane 
is given by 



Eli = <°c - <r m )KHiS, 



(A3) 



where R is the cluster radius, H± is the cluster-membrane 
height at the interface (r = R), and S is the slope of the 
cluster-membrane at the interface. This energy describes 
a line tension (~ 2irR) and vanishes for a c = a m since 
cluster and membrane contributions point in opposite ra- 
dial directions. 

We apply perturbation theory to estimate the sig- 
nificance of the surface-tension contribution. For this 
purpose, we use our previously calculated height profile, 
Eqs. [TU] and [TTJ obtained without the surface tension 
term. Using a m = a c /4 = l/c^T/mm 2 [49j and a 
physiological cluster radius R = 100 nm, we find that 
the estimated energy contribution from the surface 
tensions is much smaller than the elastic energy Eq. [21 
i.e. E^/Eei w 0.004, justifying the neglect of this term. 

Gaussian curvature 



The contribution to the elastic energy from the Gaussian 
curvature can be neglected for homogeneous membranes 
that do not change their topology (Gauss-Bonnet theo- 
rem). However, this contribution is technically non-zero 
for our cluster-membrane system. In the small deforma- 
tion approximation, this contribution is given by [39j 



+K G , 



d 2 h c d 2 h r 



d 2 h. 



dx 2 dy 2 

m W 'I'm 



dx 2 dy 2 



\ dxdy J 

(d 2 h m Y 
\ dxdy J 



d 2 r 



d 2 r, (A4) 



where K Gc ^ m ^ is the Gaussian curvature modulus for the 
cluster (membrane). In (49[, Eq. 129 shows that the 
Gaussian energy contribution has two parts. One topo- 
logical part, which is just a constant since our membrane 
contains a single receptor cluster. The second part is a 
contour integral along the cluster boundary. Based on 
Eq. 131 of [49(, we obtain 



ir(K G , c - K G , m )S 2 . 



(A5) 



For our height profile Ec 



A" 



G,c(m) 



/2 [3J 



[TU1 and [TTJ as well as 
this energy contribution is 



significantly smaller than the elastic energy in Eq. [21 i.e. 
\E^\/E e i w 0.0006, justifying the neglect of this term. 
As expected, this energy contribution vanishes when 
the elastic properties of cluster and membrane become 
identical for K GtC — K GtTn . 
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